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SECTION  I 


INTRODUCTION 

The  development  of  intense  photon  sources  for 
nuclear  weapons  effects  simulation  has  been  actively 
supported  by  the  Defense  Nuclear  Agency  for  at  least 
the  last  twenty  years.  During  this  period  it  has  become 
quite  apparent  that  pulse  power  generators  can  be  used 
successfully  to  heat  material  loads  to  the  plasma  state 
at  temperatures  and  densities  capable  of  exciting  the 
K-shell  of  most  low-Z  elements.  These  subsequently 
radiate  as  intense  X-ray  sources.  The  material  loads 
most  often  encountered  with  pulse  power  devices  are: 

(1)  an  array  of  finely  suspended  wires,  (2)  a  gas  puff, 
and  (3)  a  cylindrical  annular  foil.  These  material 
loads  are  rapidly  vaporized  by  fast  rising  mega-ampere 
currents  and,  subsequently,  converted  to  hot  dense  plasma 
that  are  copious  emitters  of  intense  X-rays.  Because  of 
the  complexity  of  the  phenomena  involving  radiation 
magnetohydrodynamics,  it  is  not  surprising  that  present 
day  understanding  of  the  plasma  behavior  leading  to  the 
emission  of  X-rays  is  still  rudimentary.  On  the  other 
hand,  considerable  progress  has  been  made  in  understanding 
and  modelling  the  radiative  properties  and  energetics  of 
low-Z  and  moderate-Z  materials.  Based  on  this  understanding 
and  on  the  empirical  information  regarding  the  temperature 
and  density  structure  across  the  radiating  plasma,  it  has 
appeared  possible  to  derive  a  relationship  which  would 
show  the  manner  in  which  the  K-shell  radiated  power  varies 
as  a  function  of  temperature,  density,  plasma  size,  and 
material  (i.e.,  Z) .  This  topic  is  the  subject  of  this 
paper . 


Before  proceeding  to  the  development  of  the  actual 
scaling  law,  it  is  reasonable  to  review  the  application 
and  importance  of  such  work.  Present  results  have 
generally  established  key  features  in  the  plasma  structure 
of  the  majority  of  high  K-shell  yield  experiments.  The 
structure  consists  of  a  set  of  compressed  hot  spots  or 
beaded  spots,  which  appear  to  emit  K-shell  radiation.  The 
spots  are  surrounded  by  an  outer  corona  which  predominantly 
radiates  at  a  cooler  temperature  in  the  L-shell.  This 
would  be  typical  of  aluminxim,  for  example.  The  temperature 
structure  suggests  a  distribution  of  ionization  states 
across  the  plasma  that  reflects  the  core/corona  pictorial 
representation  and  clearly  indicates  that  not  all  of  the 
plasma  mass  is  heated  to  the  K-shell.  At  higher  power 
levels,  it  is  also  observed  that  more  mass  is  heated  to 
the  K-shell  level,  but  that  substantially  the  same  plasma 
structure  is  maintained.  For  higher  Z  loads,  wherein 
K  =  >  L-shell  and  L  =  >M-shell,  we  still  note  that  the 
characteristic  structure  is  preserved,  i.e.,  core/corona. 
This  result  raises  an  important  question  about  the  behavior 
at  higher  power  levels;  namely,  will  higher  power  achieve 
higher  temperatures  or  will  it  simply  heat  more  of  the 
plasma  to  existing  or  slightly  higher  temperatures.  In 
the  former  case,  the  spectrvim  will  be  hardened,  although 
the  yields  may  not  increase  significantly  above  existing 
levels,  while  for  the  latter  case  the  spectrum  will  be 
preserved  but  the  yield  will  increase  substantially.  Since 
the  goal  of  the  program  is  to  both  harden  the  spectrum  and 
increase  the  yields,  it  seems  inevitable  that  a  larger 
machine  will  have  to  be  built  in  order  to  experimentally 
ascertain  this  information.  However,  from  a  theoretical 
point  of  view  it  should  be  possible  to  make  some  intelligent 


comments  regarding,  at  least,  the  scaling  of  the  yield 
in  terms  of  radiative  properties  of  the  load.  To  this 
purpose  we  derive  a  simple  scaling  law  based  on  the 
radiative  capacity  of  the  plasma  at  temperatures  capable 
of  supporting  K-shell  radiation  and  where  the  size  of 
the  emitting  region  does  not  affect,  significantly,  the 
emitting  qualities  of  the  source.  Reabsorption  through 
the  cooler  regions  has  been  neglected  for  simplicity. 


SECTION  II 


THEORY  AND  RESULTS 

Niunerical  simulations  and  spectroscopic  analysis 

of  wire  or  puff  gas  implosions  suggest  that  the  inner 

portion  of  the  compressed  plasma  is  hot  enough  to  ignite 

the  K-shell  of  all  low-Z  and  some  medium-Z  elements. 

Estimates  of  the  plasma  density  in  the  high  temperature 

19  20  3 

core  region  is  on  the  order  of  10  to  10  ions  per  cm 
depending  on  whether  the  wires  or  puffs  are  used  as  the 
load  material.  Although  these  conditions  are  not  those 
of  a  plasma  in  corona  equilibrium,  it  is  well  known  from 
the  work  of  the  Plasma  Radiation  Branch  at  the  Naval 
Research  Laboratory  that  for  the  K-shell  emitting  region 
of  the  plasma,  a  corona  equilibrium  analysis  provides 
reasonably  acceptable  results  for  the  radiated  K-shell 
energy.  Stretching  this  point  a  bit  further,  we  will 
adopt  the  assumption  of  corona  equilibrium  for  the  core 
plasma  euid  develop  a  simple  scaling  law  for  the  radiated 
power  in  terms  of  temperature,  density,  and  atomic  number. 
Opacity  effects  will  be  completely  ignored,  but  can  be 
accounted  for  in  a  post  hoc  manner. 

A  survey  of  the  literature  (see,  for  example,  the 
corona  equilibrium  calculations  by  Davis  and  Jacobs  or 
D.  Post)  reveals  that  line  radiation  is  dominant  in  hot 
plasmas  up  to  about  Z  »  40.  Beyond  Z  =  40,  bremsstrahlung 
radiation  begins  to  dominate  the  power  spectrum.  Detailed 
calculations  in  conjunction  with  the  observations  indicate 
that  about  70%  of  the  K-line  yield  emanates  from  the 
resonance  levels  of  hydrogen-and  helium-like  ionization 
stages,  on  shots  that  are  considered  good  in  that  they  are 
of  high  yields.  Using  these  empirical  facts  and  the 


assumptxon  of  corona  ionization  equilibrium,  the  peak 
K-line  radiative  power,  ,  associated  with  the 

resonance  lines  can  be  represented  as 

=  Ne  Ng  X(s  -►  q)  AE(s  -►  q)  (1 

q=l 

where  N^,  Nj  are  Njj,  respectively  and 

X(s  1)  =  X(ls^  -*•  ls2p)  ,  X(s  -►  2)  =  X(ls  -»■  2p)  are  the 
collisional  excitation  rate  coefficients  and  AE(s  ->■  q) 
is  the  excitation  energy  for  the  corresponding  transitions 
respectively.  Finally,  N  is  the  electron  density.  At 
peak  K-line  emission  the  detailed  calculations  of  Davis 
and  Jacobs  euid  Dus  ton  and  Davis  estimate  that  the 
fractional  population  densities  are 

=  0.37Nj  and  N2  “  O.SN^  .  (2 

In  addition  AE(s  1)  =0.95  AE(s  2)  where 
AE(s  -*■  2)  =  I|l  Hence 

Ae(s  -*•  2)  =  lOZ^  in  ev.  (3 


Substituting  these  estimates  into  equation  (1)  yields 


=  lOZ^  Nj^^0.35X(s  1)  +  o.5X(s  -  2) 


watts 


where  N  =  ZN-  has  been  used, 
e  1 

As  a  result  of  a  number  of  calculations  using  the 
Coulomb-Born  approximation,  the  Soviets  have  attempted  to 
parametrize  the  electron  impact  excitation  coefficient  in 
the  form  (see  Sobelman) 
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where  x  (Xo^  represents  the  upper  (lower)  state 

U  Xf 

ionization  energy,  is  the  lower  state  statistical 

weight,  6  *  AE/kT^  ,  and  A  «  20  for  Is  -*•  2p 

transitions.  This  expression,  equation  (5) ,  neglects 

exchange  effects  which  can  influence  helium-like 

transitions  by  factors  of  two  in  some  cases.  Bearing 
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this  in  mind  and  setting  x^g  “2  Ijj,  X2p  ~  J  ^  ' 

and  assuming  6  >>  C  at 

X**(ls  H.  2p)  «  4  X  10"®  -g  e"®  cmVsec.  (6) 

Z  3 

A  compilation  of  the  NRL  calculations  is  plotted  in 
Figure  1  and  suggests  the  following  dependence: 

T„  »  0.12  Z®*®®  (7) 

In  addition,  the  3 -dependence  has  been  curve  fitted  for 
Z  ^  42  with  the  result 

+  B)  e"®  a  0.012  Z^/^  (8) 

Hence,  the  total  power  can  now  be  written  as 

=  lO"^”^  Z^/^  (9) 

where  the  knowledge  that  the  resonance  lines  of  the  one- 
and  two-electron  systems  account  for  roughly  70%  of  the 
radiated  K-shell  yield  has  been  folded  into  equation  (9) . 

The  result  obtained  from  equation  (9)  is  compared  with 
the  more  detailed  NRL  models  and  as  seen  in  Figure  2  is 
in  reasonably  good  agreement,  i.e.,  to  within  about  a 
factor . slightly  less  than  two  for  4  <  Z  <  60. 

The  continuum  contribution  to  the  K-shell  yield 
comes  from  the  free-bound  and  free-free  processes.  In  the 
case  of  the  free-bound  radiative  recombination  contribution. 
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20 

Z  (nuclear) 

Rgure  1.  NRLData  Compared  to  Equation  (7) 

♦  =NRL  Data 
—s  Equation  (7) 


only  two  edges  contribute  to  the  K-line  energies  and  can 
be  represented  by 

=  1.6  X  10"^^  I  Xi  ct.  (10) 

where  x  is  the  ionization  energy,  a  is  the  rate 
coefficient,  and  N.  is  the  recombining  ion  density. 

From  a  number  of  detailed  calculations  it  can  be  shown 
that  the  following  approximations  are  reasonable, 
namely 

^bare  ~  ^He  "  ^ 

Nucleus 

If  we  take  a  =  5.2  x  10  and  substitute  the 

above  approximations  into  obtain 

=  5  X  lo"^^  Nj^  (11) 

The  bremsstrahlung  emission  is  given  by 

=  1.5  X  10“^^  Z^  Nj  J  e  d(hv) 

hv' 

where  hv'  is  the  energy  cutoff  below  which  photons  do 

not  contribute  to  the  K-shell  and  is  set  by  the  energy  of 
2 

the  Is  Is  2p  edge.  Hence 

kJ,  ^  5.3  X  10-”  z'-''  z"'"®  (12) 

Figure  3  shows  a  comparison  between  the  total  bremsstrahlung 
emission,  i.e.,  hv*  »  0  and  that  predicted  from  equation 
(12) .  The  difference  is  about  40%. 

The  total  power  radiated  at  or  above  K-line  energies 
is  shown  on  Figure  4.  The  discrepancy  between  the  two 
curves  labeled  "total  emission”  is  due  to  inclusion  of  all 


Z  (nuclear) 

Figure  4.  Radiated  Power 

---  sTotal  (all  energies) 

—  s Total  (above  K-line) 

-  •  -  s  K'iine 

-••-sBremsstrahlung  (above  K-line) 


energies  in  the  free-bound  and  free-free  estimates.  Even 
so,  the  difference  only  amounts  to  about  a  factor  of  2  at 
higher  Z  values. 

The  values  obtained  here  are  reasonable  in  comparison 
with  some  of  the  experimental  values  for  a  given  Z.  How- 
even,  this  simple  scaling  law  should  not  be  construed  as 
anything  other  than  what  it  is;  a  simple  scaling  law  for 
estimating  radiative  yields.  In  the  future,  this  model 
will  be  incorporated  into  a  simple  circuit-implosion  model. 
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ATTN 


G.  Frazier 


RAO  Associates 

AHN:  P.  Haas 
ATTN:  A.  Latter 


RAD  Associates 

ATTN:  P.  Turchl 


S-CU8ED 

ATTN:  A.  Wilson 

Science  Applications,  Inc 
ATTN:  K.  Sites 


Kaman  Sciences  Corp 
ATTN:  S.  Face 

Berkeley  Research  Assiclatas 
4  cy  AtTN:  J.  Workman 
4  cy  ATTN:  C,  Prattle 


TRW  Electronics  A  Defense  Sector 
ATTN:  Tech  Info  Ctr 
ATTN:  D.  Clement 

Kaman  Tempo 

ATTN:  DAS I AC 
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